Investigation of the spin state of Co in LaCoOa at room temperature 
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We investigate the spin state of LaCoOs using state-of-the-art photoemission spectroscopy and 
ab initio band structure calculations. The GGA-|-i7 calculations provide a good description of 
the ground state for the experimentally estimated value of electron correlation strength, U. In 
addition to the correlation effect, spin-orbit interaction is observed to play a significant role in the 
case of intermediate spin and high spin configurations. The comparison of the calculated Co 3d 
and O 2p partial density of states with the experimental valence band spectra indicates that at 
room temperature, Co has dominant intermediate spin state configuration and that the high spin 
configuration may not be significant at this temperature. The lineshape of the La 5p and O 2s core 
level spectra could be reproduced well within these ab initio calculations. 

PACS numbers: 71.20.-b, 75.20.Hr, 71.27.-fa, 79.60.Bm 



I. INTRODUCTION 

The evolution of the spin state of Co in LaCoOa with 
temperature has drawn a great deal of attention dur- 
ing last 50 years J-^i^iii^i^iii^i^iiS The ground state of this 
compound is believed to be nonmagnetic (spin S = 0) in- 
sulator. It shows two magnetic transitions at about 100 
K (a sharp transition) and at around 500 K (broad).!! 
The second transition is also accompanied by an insula- 
tor to metal transition. All these magnetic transitions 
have been attributed to the temperature induced spin 
state transition of 00^^+ ions. Initially it was believed 
that the transition at 100 K occurs due to the change 
in the spin state of Co'^+ ion from low spin (LS) state 
S* = 0) to a mixed LS and high spin (HS) 
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scopic understanding of these transitions, photoemission 
spectroscopy (PES) and x-ray absorption spectroscopy 
(XAS) have been employed extensively, since these tech- 
niques help to probe the electronic structure directlyi^i^ 
The experimental results were simulated using configura- 
tion interaction (CI) calculations for the cluster of CoOg 
octahedron. All these investigations inferred the pres- 
ence of varying mixtures of LS and HS states above 100 
K. 

Subsequently, a detailed study^ based on LDA+U 
(LDA = local density approximation and U = electron- 
electron Coulomb repulsion strength) calculations, at- 
tributed the 100 K transitions to LS to orbital ordered 
intermediate spin (IS) state {tlgel =^ S = 1) and the 
one at 500 K to orbital disordered IS state as the en- 
ergy of the latter state was found to be higher than the 
former and lower than HS state. Since then, numer- 
ous experimental and theoretical works have been car- 
ried out, which attributed the first transition to LS to IS 
state>^i^>^>^>^i^>^ This has also been demonstrated 



by PES and XAS studies,!^ where Co 2p core level, va- 
lence band and O K-edge XAS spectra of LaCoOa re- 
vealed 70% LS and 30% IS states in the temperature 
range 100 K to 300 K. No contribution of HS state was 
observed up to 300 K. 

Recently, based on GGA+J7 (GGA = generalized gra- 
dient approximation) calculations, Knizek et al. showed 
that the mixed LS-HS state (1:1) is the first excited state 
and thus attributed the first magnetic transition from LS 
to LS-HS state»^ This is in sharp contrast to the be- 
lief of IS contributions in the intermediate temperature 
range. In addition, experimental results based on XAS 
and magnetic circular dichroism^ and inelastic neutron 
scattering— also interpreted in terms of mixed LS-HS 
state in the temperature range up to 700 K considering 
cluster approximations. It is, thus, evident that the spin 
state of Co at different temperatures is still controversial. 

In this work, we investigate the spin state of Co in 
LaCoOa at room temperature using state-of-the-art pho- 
toemission spectroscopy and ab initio band structure cal- 
culations. We observe that spin-orbit coupling (SOC) 
plays an important role in determining the electronic 
structure in this system. The comparison of the high 
resolution spectra and the calculated results suggest that 
the electronic structure of LaCoOs at room temperature 
has large IS contributions. The ab initio calculations also 
provide a good representation of the experimental La bp 
and O 2s core level spectra. 



II. EXPERIMENTAL AND COMPUTATIONAL 
DETAILS 

Single crystal of LaCoOs was grown and character- 
ized as described elsewherci^ The photoemission spec- 
tra of the valence band and shallow core levels were 
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recorded at room temperature (RT) using a spectrom- 
eter equipped with monochromatic Al Ka (1486.6 eV) 
and He i (21.2 eV) sources, and a Gammadata Scienta 
analyzer, SES2002. The base pressure during the mea- 
surements was about 4x10^^^ Torr. The energy resolu- 
tion for the a;-ray photoemission was set to 0.3 eV for the 
valence band and 0.6 eV for the core level spectra. The 
resolution for the He i spectrum was fixed to 4 meV. The 
sample was cleaned in situ by scraping the sample sur- 
face using a diamond file. The cleanliness of the sample 
was ascertained by tracking the sharpness of O Is peak 
and absence of C Is peak. The Fermi level was aligned 
by recording the valence band spectrum of an Ag foil 
mounted on the same sample holder. 

The LDA-l-t/ and GGk+U {U represents the elec- 
tron correlation strength among Co 3d electrons), spin- 
polarized density of states (DOS) calculations were car- 
ried out using LMTART 6.61 (Ref. 24). For calculat- 
ing charge density, full-potential linearized Muffin-Tin 
orbital (LMTO) method working in plane wave repre- 
sentation was employed. In the calculation, we have used 
the Muffin-Tin radii of 3.515, 2.005, and 1.64 a.u. for La, 
Co and O, respectively. The charge density and effective 
potential were expanded in spherical harmonics up to / = 
6 inside the sphere and in a Fourier series in the intersti- 
tial region. The initial basis set included 6s, 5p, 5d, and 
4/ valence, and 5s semicore orbitals of La; 4s, 4p, and 
3(i valence, and 3p semicore orbitals of Co, and 2s and 
2p orbitals of O. The exchange correlation functional of 
the density functional theory was taken after H.S. Vosko 
et al. (Ref. 25) and GGA calculations are performed 
following J. P. Perdew et al. (Ref. 26). The calcula- 
tions were performed by taking LS, IS and HS configu- 
rations, which correspond to (ifgi^fgi), {t2gi^2gi^g-\) ^^"^ 
(^igT^sgX^gt) elcctronic configurations, respectively as ini- 
tial input. Self-consistency was achieved by demanding 
the convergence of the total energy to be smaller than 
10"'' Ryd/cell. Final orbital occupancies for Co t2g and 
eg states were obtained from self-consistent GGA-I-J7 cal- 
culations for different initial state configurations. (8, 8, 
6) divisions of the Brillouin zone along three directions 
for the tetrahedron integration were used to calculate the 
DOS. 



III. RESULTS AND DISCUSSIONS 

It is well known that although the electron correlation 
effects are underestimated in ab initio band structure cal- 
culations within LDA, these results are often found to 
explain most of the features in PES and XAS results in 
transition metal oxide systems j^i^i^i^Si^ However, one 
needs to consider the correlation effects to capture the 
details of the spectra and the ground state properties of 
the system. We have calculated various contributions in 
the DOS of LaCoOs using both LDA and LBA+U meth- 
ods. In order to minimize uncertainty in the calculated 
results, we have fixed the value of U to the experimen- 
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FIG. 1: (color online) O 2p (thin solid line) and Co 3d partial 
density of states having t2g (dashed line) and Eg (thick solid 
line) symmetries from (a) LDA and (b) LDA-|-(7 calculations. 

tally estimated values of 3.5 eVi^ 

In Fig. 1(a), we show the calculated LDA results for 
both the up and down spin states. There are several 
features observed in different energy positions. The fea- 
tures between -1.5 eV to -3 eV energies are dominated by 
the O 2p states with negligible contribution from other 
electronic states. These features are identified as O 2p 
non-bonding states. The features at lower energies have 
dominant contributions from O 2p electronic states with 
small but finite intensities from the Co 3d states hav- 
ing t2g and eg symmetries. Thus, the features between 
-3 to -4.5 eV can be attributed to bonding states with 
t2g symmetry and those below -4.5 eV are the bonding 
states with eg symmetry. The antibonding features ap- 
pear above -1.5 eV and have predominantly Co Zd char- 
acter. Co Zd bands with t2g symmetry appear between 
-1.5 to 0.5 eV and the eg bands between -1 to 3 eV. It 
is evident that the density of states at the Fermi level 
is large suggesting a metallic phase in the ground state 
in contrast to the insulating phase observed in various 
experiments. 

The effect of electron correlation in the density of 
states is manifested in the LDA-|-[/ results as shown in 
Fig. 1(b). Although electron correlation among O 2p 
electrons is not considered in the calculations, the en- 
ergy distribution of the O 2p partial density of states 
(PDOS) appears to be somewhat different from the LDA 
results. The O 2p contribution in the antibonding re- 
gion appears to enhance in the LDA-I-J7 results compared 
to the LDA results suggesting a spectral weight trans- 
fer from the bands having bonding character. The non- 
bonding O 2p contributions also shift toward the Fermi 
level. The changes in the Co 3c? bands are most sig- 
nificant as expected. The bonding and antibonding eg 



-6 -4-2 2 

Energy (eV) 

FIG. 2: (color online) Co 3d and O 2p partial density of states 
calculated using (a) LDA+U {U = 3.5 eV), (b) LDA+U {U 
= 4 eV) and (c) GGA+U (17 = 3.5 eV) methods. 

bands shift towards higher energies. The bonding t2g 
band appears almost at the same energies along with a 
significant increase in the spectral weight. Subsequently, 
the antibonding t2g band having primarily Co 3d charac- 
ter becomes narrower along with a reduction in spectral 
weight. It is, thus, clear that consideration of U in the 
calculations leads to a decrease in Co 3c? character and 
an increase in O 2p character of the bands in the vicinity 
of the Fermi level. Although the overlap between the t2g 
and Cg band is minimized for these parameters, no hard 
gap is observed characterizing the systems to be metallic. 

In order to investigate the value of U that creates a 
gap at the Fermi level, we compare the Co 3d and O 2p 
PDOS for U = 3.5 eV and 4 eV in Fig. 2(a) and 2(b), 
respectively. It is evident that, an increase in U to 4 eV 
generates a band gap of about 0.2 eV. The shape and en- 
ergy distribution of the DOS is almost the same in both 
the cases. The GGA is known to provide a better de- 
scription of the exchange-correlation functional as it also 
considers first order correction to the spatial distribution 
of the electronic charge density used in the LDA calcu- 
lations. In GGA, total energy of the system improves^S. 
and it is expected that the energy position of different 
bands would also improve. The calculated PDOS of Co 
3d and O 2p corresponding to GGA+U calculations are 
shown in Fig. 2(c). In contrast to the metallic phase ob- 
served in LDA results, a value of [/ = 3.5 eV creates a 
band gap of about 0.22 eV at the Fermi level in the re- 
sults corresponding to GGA calculations, which is very 
close to that observed experimentallyi^^'^ However, the 
shape and positions of Co 3d and O 2p bands is very sim- 
ilar in the LDA+U and GGA+U results. Thus, we used 
GGA+U method in the rest part of our study to discuss 
the experimental results. 
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FIG. 3: (color online) Calculated Co 3d and O 2p partial 
density of states corresponding to low spin (LS), intermediate 
spin (IS) and high spin (HS) configurations using GGA+U 
method. 

We now turn to the question of the influence of var- 
ious spin states in the electronic structure of LaCoOa. 
The PDOS obtained for different spin state configura- 
tions are shown in Fig. 3 for the same values of C/ (= 
3.5 eV). It is evident that although the LS configuration 
leads to insulating ground state, the electronic structure 
corresponding to IS and HS configuration converges to 
metallic phase. This is a well know fact for this system 
and one needs to consider additional parameters such as 
orbital ordering to achieve insulating phase.- However, 
the energy distribution of the DOS over the large bind- 
ing energy scale remains very similar. There are three 
distinctly separable features A, B and C in the DOS in 
Fig. 3. In the case of IS configuration, the t2g band be- 
comes partially filled. This effect enhances further in the 
case of HS configuration. The total intensity of the fea- 
ture A gradually decreases with the change in spin state 
configuration from LS ^ IS — > HS. In addition, the Co 
3d character of the feature A becomes highest in the IS 
configuration (80.3%) while it is about 63.8% in the LS 
state and 64.8% in the HS state. The feature B repre- 
senting the non-bonding O 2p states appears at slightly 
higher binding energies in the IS and HS states compared 
to that in the LS state. The change in the feature C is 
again substantial. The O 2p and Co 3d contributions 
are almost similar in LS and IS cases. In the case of HS 
states, the feature C has dominant O 2p character and a 
new feature D having dominant Co 3d character appears 
beyond 6 eV binding energies. 

All the above calculations are done without consid- 
ering the spin-orbit coupling (SOC) in any of the elec- 
tronic states. Recently, there is a growing realization 
that SOC plays crucial role in determining the elec- 
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FIG. 4: (color online) Co 3d and O 2p partial density of 
states corresponding to low spin (LS), intermediate spin (IS) 
and high spin (HS) configurations calculated including spin- 
orbit interactions. All these calculations are performed using 
GGA+U method. 

TABLE I: Occupancies of different Co 3d orbitals obtained 
from GGA+f/ calculations for low spin (LS), intermedi- 
ate spin (IS) and high spin (HS) configurations of Co'^"'' in 
LaCoOs. Numbers written in normal and bold face corre- 
spond to calculation without spin-orbit coupling and with 
spin-orbit coupling, respectively. 



Initial spin state 


*29T 




^2gl 




Total 


LS 


2.77 


0.58 


2.77 


0.58 


6.7 




2.77 


0.59 


2.77 


0.58 


6.71 


IS 


2.79 


1.5 


1.95 


0.43 


6.67 




2.8 


1.36 


2.03 


0.45 


6.64 


HS 


2.82 


1.9 


1.32 


0.42 


6.46 




2.83 


1.9 


1.35 


0.41 


6.49 



tronic structure in various transition metal oxides. The 
importance of SOC has been demonstrated in the sys- 
tems possessing rare-earth 4/ electrons. '^'^ In Ca2Ru04, 
it was shown that SOC plays an important role in the 
changeover of the spin and orbital anisotropy as a func- 
tion of temperaturei^ Very recent theoretical studies of 
CaaCoRhOe compound have shown that the insulating 
ground state can be achieved by considering SOC in the 
calculations.^^ In order to investigate the effect in the 
present case, we have calculated the PDOS in all the 
cases including SOC. The results are shown in Fig. 4. It 
is clear that SOC has no significant influence in the case 
of LS configuration. This is presumably due to the fact 
that all the energy bands are completely filled. The ef- 
fect is most pronounced in Fig. 4(b) corresponding to IS 
state exhibiting a splitting of the sharp feature A. 
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FIG. 5: 



(color online) Spin polarized partial density of states 
corresponding to Co t2g and Eg bands calculated (a) without 
spin-orbit coupling and (b) with spin-orbit coupling. The 
inset shows d^2_y2 and d^2 contributions in the e^j band. 



The effect of spin-orbit coupling in the occupancies of 
t2g'\ , , t2gii and Cgi are given in Table 1. Num- 
bers written in normal and bold face correspond to with- 
out SOC and with SOC, respectively. The total elec- 
tronic occupancies of Co 3c? bands corresponding to LS, 
IS and HS states are about 6.7, 6.65, and 6.47, respec- 
tively. These values are closer to those obtained from 
full- mult iplet configuration interaction calculations car- 
ried out by Saitoh et al. (Ref. 19) and somewhat smaller 
from those obtained by Korotin et al. (Ref. 9). It is 
evident that the total occupancies as well as partial oc- 
cupancies of t2g and Cg orbitals corresponding to LS and 
HS states are not affected significantly by SOC. However, 
in IS state, the occupancies of e^i orbitals decreases by 
0.14 and that of t2gi orbital increases by 0.08 due to the 
spin-orbit coupling keeping the total occupancy almost 
unchanged. 

The details of the effect of SOC in the IS state are 
shown in Fig. 5, where we show the up and down spin 
PDOS of Co 3d bands separately. In Fig. 5(a), the fea- 
ture in the vicinity of the Fermi level arises due to the 
down spin t2g states. This feature splits into two dis- 
tinctly separated features (the energy separation is about 
1.2 eV). The feature at about 1 eV binding energy repre- 
sents the doubly degenerate dxz and dyz bands and the 
non-degenerate dxy band appears above the Fermi level. 
The doubly degenerate bonding Cg^ band also splits by 
about 1 cV into dj.2_j,2 and d^^ bands as shown in the 
inset of Fig. 5(b). It is thus clear that the degeneracy 
in the t2g and eg bands are partially lifted due to spin- 
orbit interactions and the effect is most pronounced in 
the IS configuration. Since, the d bands are not com- 
pletely filled, the splitting of the t2g and eg bands will 
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TABLE II: Occupancies of 3d orbitals in intermediate spin 
configuration of Co^+ in LaCoOg. Numbers written in normal 
and bold face are corresponding to calculation without spin- 
orbit coupling and with spin-orbit coupling, respectively. 



Spin 


dyz 


dxz 


dxy 




d^2 


Up 


0.93 
0.93 


0.93 
0.93 


0.93 
0.94 


0.75 
0.86 


0.75 
0.50 


Down 


0.65 
0.92 


0.65 
0.92 


0.65 
0.19 


0.21 
0.22 


0.21 
0.23 



influence the occupancy of these bands. Here, the SOC 
induced splitting leads to lowering in energy of dyz , dzx 
and dx2_y2 bands and enhances the energy otd^y and d^^ 
bands. Thus, the occupancies of dyz , dzx and d^2_j,2 are 
expected to increase and that of d^y and dz^ to decrease. 
This is clearly manifested in the partial occupancies of 
each Co 3d orbitals given in Table 2, where the occupan- 
cies ofdyzi , dzxi and d(x^-y2'i^ increase by 0.27, 0.27 and 
0.11, respectively, and that of dxyi and dz2^ decrease by 
0.46 and 0.15, respectively. 

In order to compare the calculated results with the ex- 
perimental spectra, we show the background subtracted 
experimental valence band spectra at RT in Fig. 6. The 
spectra show three distinct features at about 1.0, 2.9 and 
5.3 eV binding energies denoted by A, B and C, respec- 
tively. The relative intensity of various features are sig- 
nificantly different in the spectra using Al Ka and He i 
sources. While the feature A is most intense in the Al 
Ka spectrum, it becomes weakest in intensity in the He i 
spectrum. Subsequently, the intensities of the features B 
and C become dominant in the He i spectrum. It is well 
known that the photoemission cross section is very sensi- 
tive to the excitation energies.'^® While the cross section 
of O 2p states is larger than that of Co 3d states in the He 
I spectrum, the relative cross section becomes opposite in 
the Al Ka energies making the cross section for Co 3d 
states significantly larger than that of O 2p states. Thus, 
the changes in intensity of various features due to the 
change in excitation energy suggests that the feature A 
is predominantly contributed by Co 3d states and the fea- 
tures B and C have large O 2p character. No significant 
feature is observed between 6 to 8 eV binding energies. 
This clearly suggests that the experimental results are 
very similar to the results obtained for IS configuration 
shown in Fig. 4(b). 

In order to bring out clarity to this comparison, we 
plot the O 2p and Co 3d PDOS corresponding to IS con- 
figuration in Fig. 6. The peak position of the feature B 
is somewhat higher (0.5 eV) in the experimental spec- 
tra compared to the calculated results. Such small shift 
in the completely filled non-bonding O 2p bands has of- 
ten been observed in the LDA results due to the under- 
estimation of the electron correlation among the O 2p 
electronsjSl The energy positions of the features A and 
C in the experimental spectra, and the large Co 3c? char- 
acter of the feature A and O 2p character of feature C 



A 

c B 




87654321 
Binding energy (eV) 



FIG. 6: (color online) Experimental valence band spectra col- 
lected at room temperature using Al Ka and He I radiations. 
The lines denote the calculated Co 3d and O 2p partial den- 
sity of states corresponding to intermediate spin configuration 
of Co including spin-orbit coupling. 

is revealed remarkably in the calculated results. Even 
the small shoulder in the binding energy range 5.5 - 7.5 
eV with relatively larger O 2p character is reproduced in 
the calculated results. These results clearly suggest that 
the electronic structure of LaCoOa at room temperature 
corresponds primarily to the intermediate spin state con- 
figuration of Co and the contribution from HS state con- 
figuration may not be so significant. These results appear 
to be different from those in the recent studies based on 
cluster approxiniationsi^i Such difference may not be sur- 
prising as in the cluster calculations, the Co 3d and O 2p 
bands are approximated to atomic energy levels. 

Finally, we investigate the applicability of ab initio 
band structure calculations in understanding the shallow 
core level spectrum. Such study will be quite interest- 
ing as core-hole potential is known to influence the final 
states of the core level photoemission. The integral back- 
ground subtracted experimental spectrum along with the 
calculated one is plotted in Fig. 7(a). Three distinct fea- 
tures visible in the spectrum. The features at 16 and 18 
eV binding energies correspond to spin-orbit split La 5p 
states as marked in the figure. The feature around 21 eV 
binding energy represent the signature of photoemission 
from O 2s levels. 

We also show the La 5p and O 2s PDOS corresponding 
to LS, IS and HS configurations in different panels of 
Fig. 7. All the three features are present in the calculated 
results corresponding to LS, IS and HS states. The shape 
of the bands in each region remains almost the same in 
each spin state. The spin-orbit splitting of the La 5p 
contributions in the PDOS is about 2 eV as also observed 
in the experimental spectra. However, the peak position 
of the bands corresponding to LS state appears at binding 
energies lower by about 1 eV in comparison to energies 
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FIG. 7: (color online) (a) Background subtracted La 5p and O 
2s core level spectra (symbols) collected at room temperature. 
The solid line represent the spectrum simulated using PDOS 
corresponding to IS state configuration. La 5p (solid line) and 
O 2s (dashed line) PDOS calculated using (b) LS, (c) IS and 
(d) HS configurations of Co. 

corresponding to IS and HS states, which are closer to 
the experimental spectra. 

In order to simulate the experimental spectrum, we 
have shifted the calculated O 2s and La 5p PDOS ob- 
tained from IS configuration by 3.1 eV and 1.2 eV, re- 
spectively and then multiplied by photoemission cross 
sections of the corresponding states<^ To account for the 
lifetime broadening and resolution broadening, O 2s and 
La 5p PDOS are convoluted with a suitable Lorentzian 
and a Gaussian, respectively. The simulated spectrum 
is shown by solid line overlapped on the experimental 
spectrum in Fig. 7(a). The reproduction of the experi- 
mental spectrum is remarkable. This suggests that one 



can indeed employ band structure calculations to deter- 
mine the core level spectrum. It is to note here that 
these calculations cannot create features appearing due 
to different screening effects expected in the final state of 
photoemission. Almost perfect reproduction of the ex- 
perimental spectrum in the present case indicates that 
final state effects may not be significant in the case of La 
5p photoemission. 



IV. CONCLUSIONS 

In summary, we have investigated the electronic struc- 
ture of LaCoOs at room temperature using various forms 
of ah initio calculations and high resolution photoemis- 
sion spectroscopy on high quality single crystal. We ob- 
serve that GGA-f [/ calculations provide a good descrip- 
tion of the ground state. Spin-orbit coupling appears to 
play a significant role in the case of IS and HS states 
of Co, which is most pronounced in the case of IS state 
of Co. The calculated Co Sd and O 2p partial density 
of states corresponding to intermediate spin state of Co 
provide the best description of the experimental valence 
band spectra at room temperature. This suggests that 
the spin state of Co at room temperature presumably has 
dominant intermediate spin configuration and that the 
contribution from high spin state is not so significant at 
this temperature. The applicability of the ah initio band 
structure calculations in understanding the shallow core 
level spectrum is also studied. The calculated exchange 
splitting of La 5p states is found to be about 2 eV, which 
is identical to the experimentally observed splitting. The 
lineshape of the calculated spectrum provides a good de- 
scription of the experimental spectrum. This suggests 
that final state effects may not be significant in the case 
of La 5p photoemission in LaCoOa. 
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